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Abstract. Heavy quarkonium decays and transitions are discussed in the framework of non- 
relativistic effective field theories. Emphasis is put on the matching procedure in the non- 
perturbative regime. Some exact results valid for the magnetic dipole couplings are discussed. 


INTRODUCTION 

In the last years the B factories, CLEO and BES have produced a large amount of new 
data for heavy quarkonium observables lU] . These data are not only interesting because 
they may signal new states or new decay or production mechanisms, but also because 
heavy quarkonium is a system that to a large extent can rigorously be studied in QCD. 
Therefore, any new understanding of it may potentially provide new insight on the non- 
perturbative dynamics of QCD. 

Heavy quarkonium, as a non-relativistic bound state, is characterized by a hierarchy 
of energy scales: m, mv and mv^, where m is the heavy-quark mass and v -C 1 the 
heavy-quark relative velocity. Whenever a system is described by a hierarchy of scales, 
observables may be calculated by expanding one scale with respect to the other. An 
effective field theory (EET) is a field theory that makes this expansion explicit at the 
Eagrangian level. To be more precise, let us call H a system described by a fundamental 
Eagrangian and suppose it characterized by 2 scales: A 3> A. The EET Eagrangian, 
=SfEFT , suitable to describe H at scales lower than A, is characterized by (1) a cut 
off A II X; (2) some degrees of freedom that exist at scales lower than /i. The 
Eagrangian =^eft is then made of all operators On that may be built from the effective 
degrees of freedom and are consistent with the symmetries of the original Eagrangian 

■^EFT = X ^ 

n 

The advantage is that, once the scale /r has been run down to A, the power counting is 
homogeneous (0„) ~ A", so that the EET is, indeed, organized as an expansion in A/A. 
Despite the EET not being renormalizable in the traditional sense, it is renormalizable 
order by order in A/A. The matching coefficients c„{A/n) encode the non-analytic 
behaviour in A. They are calculated by imposing that J^eft and describe the same 
physics at any finite order in the expansion. The procedure is called matching. Einally, 
we note that if A Aqcd, c„(A//i) may be calculated in perturbation theory, if A ~ 
Aqcd, the matching must rely on non-perturbative methods. 



Several effeetive field theories for heavy quarkonium that take full advantage of the 
non-relativistie hierarehy of seales have been developed and used over the last deeade. 
For a recent review we refer to 10] • NRQCD is the EFT that exploits the hierarchy 
A = X = mv ili. Since m ::§> Aqcd> the matching coefficients of NRQCD may 
be calculated in perturbation theory. pNRQCD is the EFT that exploits the hierarchy 
A = mv S> A = mv^ llEBi- If Aqcd ~ mv^, then the matching to pNRQCD may be still 
done in perturbation theory. We call weak coupling this regime, which may be suited 
to describe ground-state quarkonium. If Aqcd ~ then the matching to pNRQCD is 
non perturbative. We call strong coupling this regime, which may be suited to describe 
excited quarkonium states. 

The fact that EFTs may be built to describe heavy quarkonium in the strong-coupling 
regime follows from the observation that the non-relativistie hierarchy of scales survives 
also below Aqcd llli. The complication of the strong-coupling regime comes from 
the non-perturbative matching and from new scales that may arise in loops sensitive to 
Aqcd- An example is the three-momentum scale ^mAQCD discussed in |0]. Neverthe¬ 
less, many advantages remain in treating even strongly-coupled heavy quarkonium in an 
EFT framework. 

In the following we shall sketch a unified framework for the description of inclusive 
and electromagnetic decays, and radiative transitions in the framework of strongly 
coupled pNRQCD. For a treatment of inclusive and electromagnetic decay widths in 
the weak-coupling regime we refer to |0] and references therein. For a treatment of 
magnetic dipole transitions in the weak-coupling regime we refer to 


NRQCD 

NRQCD is the EFT that follows from QCD when modes of energy or momentum m 
are integrated out. The structure of the EFT Eagrangian is like Eq. ([T]) with A = m and 
A = mv ~ Aqcd- The scale mv is sometimes called soft. The degrees of freedom of the 
EET Eagrangian are quarks, antiquarks and gluons with energy and momentum lower 
than m (we neglect light quarks). 

The NRQCD Eagrangian may be written as 

"^NRQCD = -^2-/+ "^4-/+-flight, (2) 


where 


+ [y/^io^Xl. ( 3 ) 

flight = (4) 

Y is the Pauli spinor field that annihilates a heavy quark of mass m, / is the corre¬ 
sponding one that creates a heavy antiquark, iDq = ido — gT^Ag and fD = zV -f 




The term stands for the 4-fermion part of the NRQCD Lagrangian (for an explieit 
expression see, for instanee, | 0 ])- 

The eoeffieient cp is a matehing coefficient of the EFT. In Eq. Q, we have made use 
of reparameterization invariance to reduce the other matching coefficients to this one. It 
is known at two loops and may be found in El] . The 4-fermion matching coefficients 
encode the contribution of the annihilation graphs. As a consequence they develop an 
imaginary part. We refer to lll2ll for an updated list of them and for references. 

Eet us give some definitions concerning the Fock space of NRQCD. If //nrqcd is 
the Hamiltonian of NRQCD, we call |n;r, R) the eigenstates of //nrqcd. and the 
corresponding eigenvalues, r stands for the relative distance of the two heavy quarks 
and R for their centre-of-mass coordinate. Both are good quantum numbers in the static 
limit. With n we indicate a generic set of conserved quantum numbers. |n;r,R) and 
£’„(r,R; V^, V^) satisfy the system of equations: 

//nrqcd |n;r,R) = J j3rW|n;r',R')E„(r',R';V,,,V^,)5^(r'-r)53(R'-R), (5) 


(m;r,R|n;r',R') = 5„^5^(r-r')5^(R-R') 


( 6 ) 


INCLUSIVE DECAYS IN PNRQCD 

pNRQCD is the EFT that follows from NRQCD when gluons of energy or momentum 
and quarks of energy larger than mv^ and quarks of momentum larger than mv are 
integrated out. The structure of the EFT Eagrangian is like Eq. ([T]) with A = mv ~ Aqcd 
and A = mv^. The scale mv^ is sometimes called ultrasoft. In the strong-coupling regime, 
if the gluonic excitations between the two heavy quarks develop a mass gap of order 
Aqcd. then they are all integrated out from the theory. Therefore, the degrees of freedom 
of the EFT Eagrangian are only singlet quarkonium fields. The Eagrangian of pNRQCD 
is very simple: 


^pNRQCD = (^do + ^ + ^-VsjS^, (7) 

S is a non-local field, function of r, R and t, 2 0 2 in spin space and a 3 ® 3 singlet in 
colour space. The trace is taken over colour and spin indices. 

All complications go into the potential Vs, which is a non-perturbative function of r 
to be determined by a non-perturbative matching procedure. In general V 5 contains also 
an imaginary part inherited from the matching coefficients of the 4-fermion operators of 
NRQCD. The matching condition reads 

(0;r',R'|//NRQCD|0;r,R)= |^+ ^^(r'-r)53(R'-R). (8) 


The matching condition determines V 5 as a function of quantities defined in NRQCD 
(the left-hand side of Eq. (HJl). Once V 5 has been determined, one may calculate the 



solutions and Eh of the Schrodinger equation 


+ ^//(r) =Eh^h{t). (9) 

Using the optieal theorem, the inelusive decay width to light particles (l.p.) is given 
by 

= —2Im(//(0)| — =^pnrqcd|^(0))- (10) 

|//(0)) stands for a quarkonium state in the rest-frame (P = 0): 

|//(0)) = J d\ J j3i?Tr|$^,(r)Sl(r,R)}|0), (11) 

where |0) is the Fock subspace containing no heavy quarks but an arbitrary number of 
ultrasoft particles. Note that, since ^h and =SfpNRQCD have been calculated through the 
matching procedure, Eq. (fTOt provides, indeed, a practical tool to calculate the inclusive 
decay width. Explicit applications of Eq. (fTOt have been worked out in OlHIil. 


RADIATIVE TRANSITIONS IN PNRQCD 

Radiative transitions may be described in the same EET framework that we have dis¬ 
cussed so far by enlarging the gauge group to SUc{3) x t/em(l)- This means that more 
degrees of freedom have to be taken into account (photons) and more operators added 
to the EET Eagrangians. We will concentrate in the following on magnetic dipole tran¬ 
sitions iflflll . 

At the level of the NRQCD Eagrangian, magnetic transitions are accounted for by 
replacing Wq —>■ Wq — and iD ^ iD -f eegA®™ in Eq. ©, 

„em — 1 

^2-f ^2-f + ^ ■ eeQB^^xir -■ [-/D x, yr 

^em „em _ i 

+{D^, C7 ■ eenB^'^}\j/ — , y/^D'a- Xj/ 

Sm-’ 4m-’ 

em _ 1 

(tJ ■ D ■ D + D ■ eegB^^ a ■ D) t// + ■ ■ ■ + [i// ^ io^xl , (12) 

and =^iight —"flight — where the gauge fields with upperscript “em” are 

electromagnetic fields and eeQ stands for the charge of the quark of flavour Q. 

The coefficients and are new matching coefficients of the EET associated 
with the electromagnetic couplings. Again, we have made use of reparameterization 
invariance to reduce their number. All coefficients are known at least at one-loop level 






(13) 


10 . In particular, we have* 


cr = l + tf^“ = 


em 

^W1 


4 as 

‘ + 32^ + ' 
4as / 1 




m 




(14) 


Kq^ is usually identified with the anomalous magnetie moment of the heavy quark. 

At the level of pNRQCD, magnetie transitions involving ultrasoft photons are de- 
seribed by adding to the Lagrangian O the eleetromagnetie Lagrangian — 
x^/xvem/4 and a term =^ypNRQCD responsible for the coupling of the quarkonium to 
the eleetromagnetie field: 


^rpNRQCD = / d^r Tr <( S^r ■ ee^E^^S 
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+ 
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|s^c7■[rx (rxeegB"”)]}s + ^^|s^c7■eeeB‘=‘"}s + ■■■ 1. (15) 


All gauge fields are calculated in the centre-of-mass eoordinate R. The field S is under¬ 
stood as a singlet also under f/em(l) gauge transformations. 

In the eentre-of-mass of the initial quarkonium state, the power eounting goes as 
follows: ~ mv, r ~ 1 /mv, the eleetromagnetie fields assoeiated to the external photons 

go like ~ k^. The eentre-of-mass derivative V aeting on the recoiling final 

quarkonium state or emitted photon is of order ky, where ky is the energy and momentum 
of the emitted photon. 

The eoeffieients V in Eq. (fTSb are the matehing eoeffieients of pNRQCD. They 
encode high-energy eontributions to the eleetromagnetie couplings and are of the same 
nature as V 5 in Eq. Q. In the strong-eoupling regime they are determined by non- 
perturbative matehing of 5-points Green funetions involving two external quarks, two 
external antiquarks and an external photon. Eet us eonsider the matehing eondition for 


* The coefficients get also QED corrections, but they are numerically negligible. 























the 1 /m operators, it reads 


(0;r',R'|(g)(7| Jd^xyr"^a ■ eegB^’^yf + [yr ^ |0)®|0;r,R) = 


y, 


(r-VV — 
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(r ■ V^)^ I (a(1) + C7■ ( 7 |eegB®” 10) 5" (r' - r) 5^ (R' - R). (16) 
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; 3 /t>/ 


Sinee eorreetions to the state |0;r,R) involving derivatives or spins are 1/m suppressed 
(see Eq. Q), o ■ eegB^’^ effeetively behaves as the identity operator. As a consequence, 
the electromagnetic matrix element decouples in the left-hand side. From the normaliza¬ 
tion condition (|6l) it follows that 


(tB 




gB 
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_ ^em 
— Cp . 


(17) 


This is a rather remarkable result that holds to all orders in the strong-coupling constant 
and non-perturbatively. It excludes that the 1 /m magnetic coupling of the quarkonium 
field is affected by any soft contribution. A fortiori, it excludes large anomalous non- 
perturbative corrections to this coupling. Similar arguments lead to the following exact 
results at order l/m^: 


a-{rxrxB) ^2 a-B aVxE GVpxr-VE 

y^ =0, =2cf”-l, (18) 

where is the static part of the y 5 potential. The first equality follows from the fact 
that Poincare invariance protects the spin-orbit coupling nTA ll^. The second one re¬ 
markably states that to all orders in the strong-coupling constant and non-perturbatively 
the existence of an effective scalar interaction, which has been often advocated in phe¬ 
nomenological models, is excluded. The third one that those matching coefficients, like 
the one in Eq. (fTTb . get only hard contributions. 

The matching of the 1/m^ terms is more complicated. One reason is that at this order 
kinetic energy and spin-dependent corrections affect the state |0; r, R) and o ■ eegB®™ 
does not behave anymore like the identity operator. 

Once the matching has been completed, the transition width is given by: 

r d^P' d^k -T 

r«^;f'r= j (2,1)3 (2,1)3 -k-P') |g(0) ^ W'(P')y(k)] (. (19) 

where 

[//(O) ^//'(-k)7(k)] 53(P' + k) = 7(k)| - Jd^R^yp^^QCDlHiO)) . (20) 


The overline stands for the sum over the final-state polarizations and the average over 
the initial state ones. Ph = (Mh.O) stands for the four-momentum of the initial-state 
quarkonium of mass Mh- The state |/7(P)) is the state (fTTTi boosted by —V/Mh- The 
Eorentz-boost transformations may be read from | TZ , 3 . 











CONCLUSIONS 


We have diseussed in an unified framework inelusive and eleetromagnetie deeays, and 
radiative transitions of heavy quarkonium in a regime where the typieal momentum 
transfer is of order Aqcd- Noteworthy, also in this situation suitable effeetive field 
theories may be eonstrueted, systematie expansions exploited and exaet results derived. 

It seems rather unlikely that the non-perturbative matehing, onee eompleted at order 
1/m^, will support the formulas traditionally and universally used so far to deseribe 
radiative transitions at relative order and derived from phenomenologieal assumptions 
|[3,0] • This may possibly shade some light, for instanee, on the radiative transition data 
for the T system reeently eolleeted at CLEO il, whose understanding is problematie 
in many phenomenologieal models. 
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